clc;clear;

V_work = 85e3;
g_max = 65e3;

syms r

Equ = g_max == V_work/(r);

r_cond = eval(solve(Equ,r));
D_cond = 2*r_cond;

D_sheath = D_cond*exp(1)

